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We inspect tlie excitation energy spectrum of a confining string in terms of solitons in an effec- 
tive field tiieory model. The spectrum can be characterized by a spectral function, and twisting and 
bending of the string is manifested by the invariance of this function under a duality transforma- 
• tion. Both general considerations and numerical simulations reveal that asymptotically the spectral 

' function can be approximated by a simple rational form, which we propose becomes exact in the 

, Yang-Mills theory. 
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The ground state energy of a long confining string is proportional to its length L. When L decreases there will be 
corrections, and in QCD we expect for the ground state energy of the confining flux 



i?o(L) = a-L + e -I- -- + 0( — ) (1) 

> . . . . 

Here a is the string tension, e the intercept, and c = is a universal constant, the Casimir energy of zero-point 
t — fluctuations that can be computed from the Nambu-Goto action in the Gaussian approximation Unfortunately, 
a simple Nambu-Goto action becomes insufficient when we attempt to describe the excitation energy spectrum. Now 
there are additional contributions that can lead to large deformations even when the underlying displacements remain 
tiny. For example, a small torsional rotation around the axis of a long string can cause distant cross-sections to rotate 
through large angles. Similarly, if a long string is slightly bent its ends can move over considerable distances. 

In the present Letter we propose an effective field theory approach to describe the excitation energy spectrum 
of a confining string. The string appears as a soliton in the effective field theory, which enables us to account 
for contributions such as bending and twisting at the (semi)classical level. This provides us with an appropriate 
^ ' background, to systematically investigate quantum corrections. Quantitatively, our arguments are based on the 
classical action [01 



(2) 



The order parameter n is a three-component unit vector n • n = 1 and is a coupling constant that determines 
the length scale. This action is known to support closed knotted strings as solitons §| - 0- According to jsj, it 
also represents a universality class that describes SU(2) Yang-Mills theory in the long wavelength limit. Indeed, (H) 
does emerge from the Yang-Mills theory in a derivative expansion Thus it becomes natural to use its string 
solutions to inspect properties of the confining flux in the Yang-Mills theory. However, ultimately we wish to employ 
universality to propose that our results remain valid even beyond the steps that lead to (||). For this we shall combine 
duality arguments with a familiar problem in the classical theory of elasticity, the evaluation of the free energy that 
describes elastic deformations of a thin rod jrj . At least to the extent our results parallel those in the classical theory 
of elasticity, we expect their validity to reach beyond the specific details of (^). 

Classically there are exactly two intrinsically small deformations of an elastic rod that can be accompanied by 
large global displacements, twisting about its axis and bending When a straight elastic rod is twisted around its 
axis, each transverse section becomes rotated by some relative angle. We align the coordinates so that the axis of 
the rod coincides with the z-axis and denote by r the (local) twist angle, defined as the angle of rotation per unit 
length of the rod around its axis. For a constant r a rod with length L acquires a total twist -0 = tL along its length. 
Combining (|l|) with the results in we then expect that the free energy of a straight, twisted confining string has 
the form 
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-Eo + Et^^st = cj-L + e + (27r2C+^)-— + 0( — ) = n\T\-(^X + ^ + \) + 0{j^) (3) 

Here T = ^i/; counts the number of twists around the axis, C is a form-factor that characterizes the string, and 
the variable A = L/\T\ denotes length per twist. We note that (^) exhibits a A — > a/A duality, reminiscent of the 
T-duality in (super)string theories. When the number of twists T remains fixed the self-dual point A = ^/a minimizes 
the ensuing contribution which yields 

En^UQ) ~ 2V^-a\T\ + e + (4) 



In the field theory model (g) a straight twisted string is described by an axially symmetric configuration, in 
cylindrical coordinates 

sm{k(p + Tz) sn\9{r) 
n(r, (y£>, z) = I cos{kip + Tz) shi9{r) | (5) 
cos 9{r) 

where k is an integer that counts the degeneracy. We substitute (^) into (H), and by defining p = er a. dimensionless 
variable, we find for the energy of a string with length L 
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L + ^■Updpil + eD-sin'eA-'^ (6) 



We draw attention to the similarity in the functional forms of (^) and (^) which can be used to relate the parameters 
a and C to integrals over 0{p). Here we have introduced the Hopf invariant H, cxplicitcly with F — (n, dn A dn) = dA 

n = —fd^xFAA = —krL = 

so that for a straight string the Hopf invariant reduces to the number of twists TL = T that a configuration with 
length L makes around the z-axis (including the multiplicity k). 

The profile 9{p) that minimizes the energy (|^) solves the pertinent Euler-Lagrange equation. Since this equation 
depends nontrivially on the various parameters in (|^), we expect a priori that the solution also attains a nontrivial, 
nonlinear dependence on the parameters. In particular we expect that for an actual minimum energy configuration 
the functional form of in L does not stand. 

We have performed an exhaustive numerical investigation of the parameter dependence for the configuration 9{p) 
that minimizes the energy (^). The result is described by a spectral function f{k^,X). It depends on the integer fc^ 
and the dimensionless combination A = (_Le)/(27r|7i|) which measures length per Hopf invariant. Quite surprisingly, 
we find that the functional form persists: We find that to a very high degree of accuracy the energy of an actual 
straight string solution is interpolated by a simple rational spectral function 

E,nUk',X) ^ ^\n\-(aX + b+^) (7) 



where a, b, c are k^ dependent numerical coefficients. In figure 1 we plot for k — 1, together with several numerically 
computed values of the energy. The high degree of accuracy of the rational spectral function over a wide range of 
values in A and k^ leads us to propose that this simple rational form might actually be exact for the straight string 
solutions of (||) (or at least to a slight perturbation of (^ within its universality class). 

We note that (0) has a definite, manifest modular structure as it should. Namely, suppose we split a string into 
two so that the total energy remains intact. Obviously, the spectral form of the energy for both of the resulting strings 
should also coincide with (^. Since the Hopf invariant |7i| is additive and A denotes length per Hopf invariant, we 
conclude that (|^) indeed does have the requisite property under splitting and joining. 
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The definite rational form leads to additional observations. For this we first note that if A = Ac minimizes 
(0), in parallel with (||) the spectral function in (|^) exhibits a A ^ A^/A duality. In order to interpret this duality 
we consider a string with Dirichlet-type boundary conditions that clamp the ends which keeps Ti intact. According 
to the string has a tendency to adjust its length A towards A = Ac i.e. towards the self-dual point of the spectral 
function that minimizes the energy: When A < Ac the value of A has a tendency to increase, but when A > Ac there 
is a tendency for A to decrease. 

We recall the relation between Hopf invariant 7i, twist T and writhe W |^ 

n ^ T + W (8) 



From this we conclude that for A ^ Ac the string has a tendency to supercoil: When A < Ac it becomes energetically 
favourable for the string to dilate. This tends to decrease twist T, hence there is a tendency to increase writhe W. 
Similarly, when A > Ac it becomes favourable to contract which can lead to an increase in twist at the expense of 
writhe. As a consequence the A ^ Ac/A duality of the rational spectral function becomes related to supercoiling. It 
can map a supercoiled configuration with an excess amount of twist to a (dual) supercoiled configuration with an 
excess amount of writhe but with an equal amount of energy: Since the Hopf invariant Ti. does not change when we 
exchange T W, we have a twist- writhe duality mapping that exchanges T ^ W and sends A — > A^/A but leaves the 
the spectral function intact. The energy (M) is invariant under this twist-writhe duality transformation, it determines 
a symmetry in our effective field theory M. 

The previous arguments are based on the inspection of a straight string. But a configuration with a nontrivial 
writhe can not be straight, it must be bent. For this we recall from classical theory of elasticity that the free energy for 
bending of a thin elastic rod, uniformly bent over a length i, has the form E = kL/ R~~^ with R the (mean) radius of 
curvature and k a form factor that characterizes the rod . Suppose we consider a straight confining string with Hopf 
invariant Ti. and a constant rate of twist along the string. We bend the string slightly in a uniform and non-planar 
manner, with a constant radius of curvature R. This leads to a decrease in twist and, since Ti. remains intact, to a 
nontrivial writhe which we describe by a constant rate of increment w along the string. For small to the length L{uj) 
of the string scales in proportion to uj according to L{uj) = Lq + c ■ uiR, where Lq is the (planar) component of the 
length that does not contribute to the writhe and c is a numerical constant; for simplicity we set c = 1. When we 
combine the free energy for bending with the ensuing free energy (||) for stretching and twisting, we conclude that to 
a leading order the free energy of our string acquires the following lo dependence, 

E{.) = aL(u.) + • + ' ^ (9) 



with a and (3 some parameters. 

We first consider the minimum value of the free energy (|^) when wc minimize it w.r.t. the rate of increment uj by 
keeping the Hopf invariant Ti fixed. We get 

E^,^ ^2{p-a)-^-(\n\ + ^\ (10) 



R \ R . 

where we have defined = + fP + aR^. We observe that (|l^) has the same functional form as (^. This suggests 
in particular, that the underlying dual structure persists. For this we subject the free energy to a minimization of w, 
with the condition that both the Hopf invariant Ti and the total length of the string L = La + loR remain fixed. We 
find 

^- ' ' M)^-^) (11) 



V \Ti\ ^ ' ' S 



where S*^ = a'^{L — Lq)'^ + 0^1/^. The dual structure is now manifest. We also note that both ( po|) and ( pT|) are 
consistent with the expected modular structure under splitting and joining. 

In the field theory model (^ a stable, bent and twisted confining string is (locally) described by a closed knotted 
soliton 1^, The radius of curvature of an actual soliton in general fails to be constant, but since the deviations 
from an average R appear to remain small B we can expect the main features of (0) to persist. The energy of a closed 
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knotted soliton should then admit an asymptotic expansion in terms of modular invariant quantities such as Hopf 
invariant and length/average value of the radius of curvature. In particular, we expect that the leading contribution 
in this expansion has the manifestly twist-writhe dual functional form (0) , (|ll|) of the spectral representation, with the 
various parameters now integrals of n evaluated for the soliton. Since the soliton minimizes the energy it corresponds 
to the self-dual minimum energy point of the spectral function. According to , ( pO| ) this means that to the leading 
order in an asymptotic expansion the energy spectrum of a closed knotted soliton admits the now-familiar functional 
form 

EsoHton = 7- + e) (12) 

with some parameters 7 and s In particular, since \H\ — 1, 2, 3, ... is an integer this energy spectrum is equipar- 
titioned. En = 7 • {n + e). 

Due to the high complexity of the Euler-Lagrange equations of (||) it becomes difficult to actually check the accuracy 
of the asymptotic (|l^), or to inspect the underlying dual structure. Numerical simulations are very demanding [Q, and 
at the moment we lack computational resources to perform exhaustive simulations. Thus we rely on the numerical 
results that have been recently published in There, the energy of various knotted solitons is evaluated for 
Qh = 1, ...,8. In figure 2 we compare the data in Q to the predictions of When we account for the numerical 
uncertainties in the data as described in [Q, we conclude that the agreement appears quite satisfactory. In particular, 
this suggest that for (Q) the dual twist- writhe structure is inherent. 

In conclusion, we have inspected properties of a confining string by describing it as a soliton in an effective field 
theory. In particular, we have numerically studied the excitation energy spectrum of a straight string. We find that 
it can be described by a spectral function that admits a simple rational form, manifestly invariant under a duality 
transformation that relates twist and writhe. We have verified that the presently available 3D data is also consistent 
with the simple rational form of the spectral function. This suggests that (||) is at least in the same universality class 
with a Lagrangian for which a manifestly dual, simple rational form of the spectral function is exact. It becomes 
natural to expect that our rational realization of the twist-writhe duality is distinct to an effective Lagrangian that 
accurately describes the properties of the confining flux in the Yang- Mills theory. 
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Figure Caption 



Figure 1: A comparison of the rational function with the energy values of the line solitons in (^. The 
interpolation is in the sense of least square, which yields a w 4.15..., h w 5.26... and c k, 1.42.... The agreement is 
consistent with finite lattice size errors. 

Figure 2: Comparison of the linear trajectory ( p^ ) (see also ||l^) to numerics in Q|. Line/data with + denotes 
the toroidal 7i = 1, 8 solutions with slope/intercept 364/116 and line/data with o the minimal energy 7i = 1, 8 
solutions with slope/intercept 277/306. (We note that considered only interpolation with a vanishing intercept.) 
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